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The topological invariant of a topological insulator (or superconductor) is given by the number of 
symmetry-protected edge states present at the Fermi level. Despite this fact, established expressions 
for the topological invariant require knowledge of all states below the Fermi energy. Here, we propose 
a way to calculate the topological invariant employing solely its scattering matrix at the Fermi level 
without knowledge of the full spectrum. Since the approach based on scattering matrices requires 
much less information than the Hamiltonian-based approaches (surface versus bulk), it is numerically 
more efficient. In particular, is better-suited for studying disordered systems. Moreover, it directly 
connects the topological invariant to transport properties potentially providing a new way to probe 
topological phases. 

PACS numbers: 72.20.Dp, 73.43.-f, 74.20.Rp 



I. INTRODUCTION 

Given a Hamiltonian H(k) of a band insulator or a su- 
perconductor and its symmetries as a function of the mo- 
mentum fc in (i-spatial dimensions, a topological invari- 
ant Q(H) can be defined. It counts the number of surface 
states insensitive to disorder which are present at an in- 
terface between the system and the vacuum. In each spa- 
tial dimension exactly 5 out of 10 Altland-Zirnbauer sym- 
metry classes (distinguished by time-reversal symmetry 
T, particle- hole symmetry V , and chiral/sub-lattice sym- 
metry C) 1 allow for a nontrivial topological invariant. 2,3 

The evaluation of the topological invariant conven- 
tionally involves an integral over a <i-dimcnsional Bril- 
louin zone of some function of the Hamiltonian. Re- 
cently, various approximations to the topological invari- 
ant have been developed which require only the knowl- 
edge of eigenvalues and eigenvectors of the Hamiltonian 
at one point in momentum space (rather than in the en- 
tire Brillouin zone). 4-6 

Despite the fact that these approximations are more 
efficient, we argue that they do not use one important 
property of a topological invariant. By definition, the 
topological invariant describes the properties of the sys- 
tem at the Fermi level, namely the number of edge states. 
This observation suggests that it should be possible to 
calculate the topological invariant without knowing the 
full spectrum of the Hamiltonian, but rather calculating 
only properties of the system at its Fermi energy. For 
one-dimensional (ID) systems, this was demonstrated in 
Ref. 7. Here, we show that for any dimensionality the 
topological invariant can be obtained from the scattering 
matrix of the system at the Fermi level. 

Our results offer two benefits. Firstly, since the scat- 
tering matrix contains less degrees of freedom than the 
Hamiltonian, the computation of the topological invari- 
ant is much more efficient. Secondly, the scattering ma- 
trix relates the topological invariant to transport prop- 
erties, suggesting ways to probe the topological phase by 
electrical or thermal conduction measurements. 8 ' 9 

The approach is based on dimensional reduction: We 



relate the scattering matrix in d dimensions to a Hamilto- 
nian in d — 1 dimensions. Our scheme of dimensional re- 
duction does not preserve the symmetry, unlike the field 
theory based scheme of Ref. 10. Instead our dimensional 
reduction preserves the topological invariant, similarly to 
the dimensional reduction of clean Dirac-likc Hamiltoni- 
ans of Ref. 11. 

In the remainder of the introduction we first illustrate 
our approach by revisiting the familiar example of the 
integer quantum Hall effect. Subsequently, we present a 
brief outline of the paper. 



Dimensional reduction in the quantum Hall 
effect 



A 2D system exhibiting the integer quantum Hall ef- 
fect is a topological insulator in the symmetry class A 
(all symmetries broken). It is characterized by a quan- 
tized transverse conductance o~ xy = ngo with n £ Z 
and go = e 2 /h. The quantum number n is a topo- 
logical invariant (the so-called Chern number) of the 
Hamiltonian. 12 It equals the number of protected chiral 
edge states at the Fermi level, each of which contributes 
e 2 /h to the transverse conductance. 13,14 

Charge pumping provides an alternative way to relate 
the topological invariant to a quantized transport prop- 
erty: inserting a flux quantum inside a quantum Hall 
sample rolled-up to a cylinder adiabatically pumps n elec- 
trons across the sample. 15 There exists a scattering ma- 
trix formulation of charge pumping, 16,17 which allows to 
express pumped charge per cycle (in units of e), 

Q = J- f ^^-logdetr(^), (1.1) 
liri J dip 

through the flux dependence of the reflection block r(ip) 
of the scattering matrix of one lead. 18 Here ip denotes the 
dimensionless flux $ = fup/e and the system is assumed 
to be insulating such that the reflection matrix r(tp) is 
unitary. Equation (1.1) is nothing but the winding num- 
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ber of det r(ip) when tp is varied from to 27r, which is a 
topological invariant. 

The winding number occurs in a different context in 
the theory of topological insulators. The topological 
invariant of a one-dimensional Hamiltonian H(k) with 
chiral/sub- lattice symmetry 



(1.2) 



19,20 



is expressed via the winding number given by 

1 r 27T d 
Q(H) = — dk — logdetft(fc). (1.3) 

Z7TI Jq are 

Here momentum k is measured in units of ft/a, with a 
the lattice constant. We see that upon the identification 
h = r and k = ip we are able to express the topological 
invariant in a 2D system without any symmetries as the 
topological invariant of an effective Hamiltonian in ID 
with chiral symmetry. We will show that a similar di- 
mensional reduction applies to all topological invariants 
in all dimensions. 



B. Outline of the paper 

As a prerequisite for the dimensional reduction, we 
have to open up the system to obtain a scattering ma- 
trix from a given Hamiltonian. Section II explains how 
this can be done. This section may be skipped on first 
reading. The dimensional reduction proceeds along the 
following lines: First we form out of a scattering ma- 
trix S a reflection block r(k) from one surface of the 
system, when all the dimensions except one are closed 
by twisted periodic boundary conditions. Then, the ef- 
fective Hamiltonian H d -i(k) in one dimension lower is 
defined according to the simple rule 



Hd-\{k) = r(k), with chiral symmetry, 



(1.4a) 



H d -i(k) = 



_ ( r{k) 



r\k) 



, without chiral symmetry. 

(1.4b) 



In Sec. Ill we show how to evaluate r(k) given the scat- 
tering matrix S of the initial system and prove that the 
reduced Hamiltonian Hd—i has the same topological in- 
variant as the original iJ, i.e. Q(Hd-i) = Q{H). 

After the general proof we turn to the particular ways 
to evaluate the topological invariant in 1-3 dimensions 
in Sec. IV. In ID we show that our expressions coincide 
with the ones derived in Rcf. 7 in a different way, with- 
out using dimensional reduction. For 2D we formulate 
the evaluation of the topological invariant as a general- 
ized eigenvalue problem. For 3D topological insulators in 
class All the topological invariant reduces to a product 
of 2D invariants, while the other symmetry classes re- 
quire usage of a Bott index. 6 We also mention how weak 
topological invariants fit into our approach. 



We consider the numerical efficiency of our method and 
show examples of its application in Sec. V. We also com- 
pare the finite size effects of different approximations to 
the topological invariant, and introduce the 'fingerprint' 
of phase transitions between different topological phases 
in 2D. Finally, we conclude in Sec. VI. 



II. SCATTERING MATRIX FROM A 
HAMILTONIAN 

This section contains the necessary preliminaries: the 
definition of scattering matrix and a proof that the shape 
of the Fermi surface can be calculated from the scattering 
matrix. 

While the formulas in this section are needed for the 
actual implementation of our method of dimensional re- 
duction, the method itself can be understood without 
them. This section can thus be skipped at first reading. 

Any Hamiltonian H(k) of a translationally invariant 
system with a finite range hopping can be brought to the 
tight-binding form by choosing a sufficiently large unit 
cell 



H(k) = H 



d 

E 

i=0 



tie 



iki 
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(2.1) 



Here k is a d-dimensional vector of Bloch momenta, H is 
the on-site Hamiltonian, and ti are the hoppings in posi- 
tive i-direction. We start our consideration from opening 
the system and attaching 2d fictitious leads to it. First 
we attach d sites to the original system without on-site 
Hamiltonian, and connect them with hoppings ti to the 
system. The Hamiltonian of this 'unfolded' system be- 
comes 



H 



H it 
t 



t — (ti,t 2 , . . . tdY 



(2.2) 
(2.3) 

In the next step we attach the fictitious leads to the 
unfolded system, as illustrated in Fig. 1 for the case of 
two dimensions. The hopping to the leads in positive i- 
direction is chosen to be equal to +1, and in the negative 
i-direction to be equal to —1. 

We are now ready to construct the scattering matrix 
of the open system by using the Mahaux-Wcidcnmullcr 
formula 21 (see also Appendix A) 

S = l + 2niW t (H -inlVW^^W. (2.4) 

The coupling W between the lead and the system is equal 
to Wy/p, with w the hopping from the lead to the system, 
and p the density of states in the lead. We choose p = 
1/wtt, such that 



W = 



/I 
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(2.5) 
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FIG. 1: Sketch of the tight binding model used to attach 
leads in order to open-up the Hamiltonian H(k) of (2.1). In 
2D we introduce four leads shown as circles labeled by 1, 1, 2, 
and 2. The on-site terms (boxes) are connected by hoppings 
(lines). The additional trivial hoppings 1 and —1 are intro- 
duced such that the lead properties drop out when twisted 
periodic boundary conditions are applied. For the Mahaux- 
Weidenmuller formula (2.4) , the three nodes form the on-site 
Hamiltonian H which is then connected via the trivial hop- 
pings to ideal leads. 



here, we have set w = 1 for convenience. The values of 
hopping and the lead density of states are chosen such 
that in the process of rolling-up, the fictitious leads drop 
out. 

The scattering matrix (2.4) relates the incoming states 
in the leads to the outgoing ones: 













0i 


"02 




02 


"02 


= s 


■02 



0, 

w 



4>d 

w 



(2.6) 



To prove that the scattering matrix contains all of the 
information about the Fermi level at energy Ep = 0, 
we impose twisted periodic boundary conditions on the 
scattering states: 
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(2.7) 
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with the twist matrix Z(k) given by 

/ e lfel ••• \ 



Z(k) 



e~ ikl 







e 
-ikd o J 



(2.8) 



5 = 1 + 2iW* 
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We show that Eqs. (2.6) and (2.7) have a solution for 
a given k if and only if the equation H(k)ip = has 
a nontrivial solution. The condition for the nontrivial 
solution of Eqs. (2.6) and (2.7) to exist is 

det[5 - Z{k)] = 0. (2.9) 

Performing block-wise inversion of H — iirWW^ yields 

J = (H Q - id - ith)- 1 . (2.11) 

We simplify this expression further by noting that 

S = lz + 2i lz uUu, U=(l-itf)W, (2.12) 

with 7 Z the third Pauli matrix in the direction space. We 
now write 

dct[S* - Z{k)] = dct [1 + lz Z(k) + 2iUljU] (2.13) 

= dct Jdet[l +7 2 Z(fe)] 

x det (J- 1 + 2iU[l + lz Z{k)]U^) 

= dct Jdct[l +j z Z(k)]detH(k). 

Since both J and 1 + r y z Z(k) are nonsingular, the last 
identity means that det [S — Z{k)\ and det H(k) can only 
be zero simultaneously, which is what we set out to prove. 

This proof shows that the Fermi surfaces as defined 
by the original Hamiltonian and the scattering matrix 
are identical. This is the reason why it is at all possible 
to determine the topological invariant using solely the 
scattering matrix S. Even though the scattering matrix 
only describes scattering at the Fermi level, it contains 
information about the complete Brillouin zone, and thus 
cannot be obtained from a long wavelength or low energy 
expansion of the Hamiltonian, but requires the complete 
Hamiltonian. Note however that the scattering matrix 
at a single energy contains less information about the 
system than the Hamiltonian: in order to determine the 
Hamiltonian from the scattering matrix, the inverse scat- 
tering problem has to be solved which requires knowledge 
of the scattering matrix at all the energies. 

The size of the scattering matrix (2.4) is 2d-tim.es 
larger than the size of Hamiltonian. However, if the 
Hamiltonian is local on a large d-dimcnsional lattice with 
size L d , the hoppings ti are very sparse. This allows to 
efficiently eliminate all of the modes except the ones that 
are coupled to the hoppings. The resulting scattering ma- 
trix is of size 2dL d ~ 1 , and accordingly for large systems 
it is a dense matrix of much smaller dimensions than the 
Hamiltonian. 
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III. DIMENSIONAL REDUCTION 

The aim of this section is to provide a route to the 
topological classification of scattering matrices by elim- 
ination of one spatial dimensions. This approach of di- 
mensional reduction is inspired by the transport prop- 
erties of topological systems. When applied to ID sys- 
tems it reproduces the results of Rcf. 7, and in quantum 
Hall systems it reproduces the relation between adiabatic 
pumping and the Chern number of Refs. 15,18. 

We begin from substituting the first 2(d — 1) equa- 
tions from (2.7) into (2.6). This is equivalent to applying 
twisted periodic boundary conditions to all of the dimen- 
sions except the last one, which is left open. Then we 
study the reflection from the d-direction back onto itself. 
The reflection is given by 

Vtyout = r(k)tp dM , (3.1) 
r(k) =D-C[A-Z d _ 1 {k)]- 1 B 1 (3.2) 

with Z<i-i given by Eq. (2.8) in d — 1 dimensions. The 
matrices A, B, C, and D are sub-blocks of S given by 
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C = (S d ,i ■ ■ ■ S^dpx) , 

To study topological properties of r(k) we construct 
an effective Hamiltonian Hd-i(k) which has band gap 
closings whenever r(k) has zero eigenvalues. In classes 
possessing chiral symmetry one may choose a basis such 
that r(k) — r T (fe). If chiral symmetry is absent, there is 
no Hcrmiticity condition on r, so we double the degrees 
of freedom to construct a single Hermitian matrix out of 
a complex one. The effective Hamiltonian is then given 
by 



Hd—i(k) = r(fe), with chiral symmetry, 

r(k) 
^(k) 



H d -i(k) = 



(3.4a) 

without chiral symmetry. 

(3.4b) 



It is straightforward to verify that in both cases the 
Hamiltonian Hd-i(k) has band gap closings simultane- 
ously with the appearance of vanishing eigenvalues of 
r(fc). 

If r(fe) has chiral symmetry, Hd-i(k) does not have it. 
On the other hand, if r(k) has no chiral symmetry, then 



H d -i{k) = -T z H d ^x(k)r z 



(3.5) 



with t z the third Pauli matrix in the space of the dou- 
bled degrees of freedom. This means that in that case 
Hd-i(k) acquires chiral symmetry. 

The way in which the dimensional reduction changes 
the symmetry class is summarized in Fig. 2. The trans- 
formation of symmetries of r(k) into symmetries of 
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TABLE I: (Color online) Topological classification of the ten 
symmetry classes in different dimensions. Combinations of 
symmetry class and dimensionality which support non-trivial 
topological invariants are indicated by the type of the topo- 
logical invariant (Z or Z 2 ). Classes which support only triv- 
ial insulators are denoted by The arrows indicate the 
change of symmetry class upon dimensional reduction as dis- 
cussed in the main text. The topmost symmetries A and AIII 
(which do not have any anti-unitary symmetries) transform 
into each other, whereas the remaining 8 classes (with anti- 
unitary symmetries) exchange cyclically. The dimensional re- 
duction changes the symmetry class, but preserves the topo- 
logical invariant ('-', Z, or Z 2 ). 



Hd-i(k) is straightforward in all of the cases, except the 
time-reversal symmetry in symmetry classes All and AI. 
There we have r(k) = ±r T (— k), and hence 



Hd-i(k) 



r(fc) 
rt(fc) 

±r(-k) T 
±[r T (-fc)]t 



±T x H* d _ l {-k)T x . (3.6) 



The details of the symmetry properties of r and H , as 
well as the relations between these symmetries are given 
in App. A. 

The way the symmetry class of the d-dimcnsional 
Hamiltonian transforms into the symmetry class of 
Hd-i(k) expresses the Bott periodicity of the topolog- 
ical classification of symmetry classes. 2 Namely, symme- 
try classes A and AIII transform into each other, and the 
other 8 classes with anti-unitary symmetries are shifted 
by one, as shown in Table I. This reproduces the nat- 
ural succession of symmetry classes that appears in the 
context of symmetry breaking 22 (see also Appendix A). 
The combined effects of the change in dimensionality and 
in symmetry class is that the Hamiltonians H(k) and 
Hd-i(k) have the same topological classification. 

We now turn to prove that for localized systems topo- 
logical invariants Q(H) and Q(Hd-i) arc identical. This 
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C ::^(k) = -r z H{k)r z 



AIII 



7> : r(k) = -r*(-k) 



t 



CI V:\r(k) = -a y r*(-k)a y 



AI 



BDI 



T : r(k) = -a y r T (-k)a y 



T 



T:[H(k) = T x H'(-k)r x r 
V : H(k)'=-r y H : (^k)r y 



•(k) = r T (-k) 



T: H{k) = r z H*{-k)r z 
V : \H\k)'=-'H*(-k)\ 



T :\H(k) = H*(-k)\ 



V : r(k) = r*(-k) 



T : r(k) = r T (-k) 



FIG. 2: Symmetry properties of r(fe) and H(k) in the ten symmetry classes. Time-reversal symmetry is denoted by T, particle- 
hole symmetry by V. The signs at the top and left of the table denote either the absence (x) of a corresponding symmetry, 
or the value of the squared symmetry operator. The entries of the table with a gray background have an additional chiral 
symmetry C, which always has the form shown in the AIII entry of the table. In particular, we always chose a basis such that 
r(fc) = j-^(fe) in the chiral symmetry classes. The way symmetry classes transform under our definition of Hd-i, cf. (3.4), is 
denoted by the arrows; the double arrow implies a doubling of degrees of freedom as in Eq. (3.4b). Going along an arrow, 
the symmetry of the reflection block r(k) (marked by a solid box) transforms into the symmetry of the reduced Hamiltonian 
(marked by a dashed box). In the chiral classes, there is an additional symmetry (not marked by a box) which can be obtained 
from the other by combining it with the chiral symmetry, H(k) = —r z H(k)r z and r(fe) = J* (fe), respectively. 



correspondence was proven in ID in Ref. 7, so here we 
accomplish the proof in higher dimensions. 

First of all, we observe that a topologically trivial 
Hamiltonian can be deformed into a bunch of completely 
decoupled localized orbitals without closing its gap. In a 
sufficiently large system, this also means that the gap of 
Hd-i(k) does not close during this process. For a system 
of decoupled orbitals, r(fe) and accordingly Hd-i(k) are 
momentum-independent (and hence Hd—i{k) is topologi- 
cally trivial) . This means that a sufficiently large system 
with trivial H(k) maps onto a trivial Hd-i(k) under the 
scheme of dimensional reduction outlined above. 

Let us now consider an interface between two systems 
with different bulk Hamiltonians H and H' , shown in 
Fig. 3. If the Hamiltonians Hd-\ and H' d _ 1 constructed 
out of reflection blocks of the two systems have differ- 
ent topological invariants, a topologically protected zero 
energy edge state in d — 1 dimensions must appear at 
the interface between them. Recalling that a zero energy 
edge state in d — 1 dimension corresponds to a perfectly 
transmitting mode of the original <i-dimcnsional system, 
we conclude that H and H' have different topological 



invariants. 

Conversely, if H and H' have different topological in- 
variants, there exists a transmitting mode at the interface 
between two parts of the system, which appears irrespec- 
tive of system size and microscopic details of the inter- 
face. This means that it is not possible to construct an 
interface between H^-i and H' d _ 1 which would be com- 
pletely gapped. 

Finally, the edge states in d — 1 dimension have to 
have the same group properties as the surface states 
in d dimensions, leading us to the conclusion that 
Q{H) = Q(Hd-i), as we set out to prove. The topology- 
preserving property of our dimensional reduction proce- 
dure is the same as that of the mapping from a general 
d-dimcnsional Hamiltonian to a d+ 1-dimcnsional Hamil- 
tonian presented in Ref. 23. 

At this point one might wonder why we apply the 
dimensional reduction only once. Indeed, the reduced 
Hamiltonian Hd-i can be straightforwardly approxi- 
mated by a tight-binding Hamiltonian on a d — 1 dimen- 
sional lattice using a Fourier transform. This allows to 
repeat the procedure of dimensional reduction until we 
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FIG. 3: A system in d dimensions consisting out of two parts 
with different Hamiltonians H and H . Reflection blocks of 
the scattering matrix r and r' are used to define the lower 
dimensional Hamiltonians Ha-i and H' d _ 1 . We prove the 
correspondence between topological invariants in d and d—1 
dimensions using the relation between the surface state at 
the interface between H and H' and the edge state at the 
interface between Hd-i and H' d _ 1 . 



arrive at a zero dimensional Hamiltonian. We stop at the 
first dimensional reduction for practical purposes, since 
the advantage of considering only Fermi level properties 
is achieved already at the first step. 



IV. RESULTS FOR ONE-THREE DIMENSIONS 

A. Topological invariant in ID 

We begin by verifying that we recover the ID results of 
Ref. 7, where the topological invariant was related to the 
scattering matrix without going through the procedure 
of dimensional reduction. Dimensional reduction in this 
case brings us to a zero-dimensional Hamiltonian. The 
topological invariant of a zero-dimensional Hamiltonian 
without symmetry between positive and negative ener- 
gies (symmetry classes A, AI, and All) is given just by 
the number of states below the Fermi level. In class All 
Kramers' degeneracy makes this number always even. In 
addition, in OD there exist two Z2 topological insulators 
in symmetry classes D and BDI. The topological num- 
ber is in that case the ground state fermion parity, or 
the Pfaffian of the Hamiltonian in the basis where it is 
antisymmetric. To summarize, 

Q(H) = u(H), for A, AI, and All, (4.1a) 
Q(H) = Pf iH, for D and BDI, (4.1b) 

where v(A) denotes the number of negative eigenvalues of 
the Hcrmitian matrix A. Substituting H from Eqs. (3.4) 



2 = v(r), 
Q = Pf ir, 

Q = Pf 



for AIII, BDI, and CII (4.2a) 
for Dili, (4.2b) 



ir 
-ir T 



det r, 



for D. 



(4.2c) 



We confirm that the Eqs. (4.2) arc in agreement with 
Ref. 7. 



B. Topological invariant in 2D 

Starting from 2D, the dimensional reduction brings us 
to a ID Hamiltonian. In this subsection we first review 
the known expressions for the topological invariants of 
ID Hamiltonians, and then describe how to efficiently 
evaluate it for the effective Hamiltonian (3.4). The Z 
topological insulators in ID (classes AIII, BDI, and CII) 
are characterized by a winding number 19 ' 20 



Q(H) = ^- [ log det /*(&), 

2m J ak 



(4.3) 



for AIII, BDI, and CII. (4.4) 

The topological invariant for the Hamiltonian in class D 
is given by Kitaev's formula 24 



Q{H) = sign 



Pf H{0) 



Pf H{n) 



for D. 



(4.5) 



Finally, in class Dili the expression for the topological 
invariant was derived in Ref. 25: 



Q(H) 



■ exp 



Pf[U T h(0)\ 

¥f[U T h{-K)] Vdct h(0) 
Pf [U r h(0)} ^det h{n) ' 



dk — log det h(k) 
ak 



for Dili, 



(4.6) 



where the square root is defined through analytic con- 
tinuation over the first half of the Brillouin zone, h is 
defined by Eq. (4.3), and Uj- is the unitary part of the 
time reversal operator T = Uq-JC. 

Substituting Eq. (3.4) into the expressions for topolog- 
ical charge we get 

1 f 2 ^ d 

Q = — / dk— logdetr(fc), for A, C, D (4.7a) 
2m' J Q dk 



Q = 
Q = 



Pf[U T r(ir)] Vdetr(O) 
Pf[Urr(0)] ydetr^)' 
Pf r(0)" 



sign 



Pf r{n) 



for All, (4.7b) 
for DHL (4.7c) 
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In order to efficiently evaluate the integral given in 
Eq. (4.4), and the analytic continuation in Eq. (4.7b) 
using Eq. (3.2), we define a new variable z = e lk . Then 
we perform an analytic continuation of det r{z) to the 
complex plane from the unit circle \z\ = 1. To find zeros 
and poles of det r(z) we use 



det r(z) = det 



A-Zx(k) B 



C 



D 



dctL4-Zi(fc)] 



where 



(4.8) 



Zi(k) 



ik 



e 
-ik o 



Equation (4.8) follows from Eq. (3.2) and the determi- 
nant identity 



det(D - CM~ l B) = det 



M B 
C D 



det M. (4.9) 



Since both the numerator and the denominator of 
Eq. (4.8) are finite at any finite value of z, the roots of 
the numerator z n are the zeros of detr(z), and the roots 
of the denominator w„ are the poles. In App. B we show 
that due to unitarity of the scattering matrix, the poles 
of detr(z) never cross the unit circle. By multiplying 
the second column of the numerator of Eq. (4.8) by z we 
bring the problem of finding roots z n of this numerator 
to the generalized eigenvalue problem, 



2.2 



which can be efficiently evaluated. The roots w„ of the 
denominator can also be found by solving the generalized 
eigenvalue problem, 



£1,1 


-1 S 


Si.! 


s 


$2,1 


s 



\ 








Sin 




/ v° 


~ $2,1 





S1.1 —1 1 
5n )^ = w - 



-s hl 

1 Si,i 



Ipn- (4.11) 



Since the poles of det r(z) never cross the unit circle, 
in classes A, C, and D the topological invariant is given 
by 

Q = #{z n : \z n \ < 1} - iVi, for A, C, and D, (4.12) 

i.e., the number of z„'s inside the unit circle minus the 
number of modes N\ in the direction 1. In class All 
(quantum spin Hall insulator) the topological invariant 
is given by 



Q = 



n, 



Pf U r r{n) 
Pf U T r{0) 



for All, (4.13) 



with the branch cut of the square root along the negative 
real axis. Note that the linear fractional transformation 
z 1 y i (1 + z) I ( 1 — z) maps the upper half of the unit circle 



onto the negative real axis. In symmetry class Dili the 
evaluation of the topological invariant is most straight- 
forward, and yields 

"Pf r(0)' 



Q = sign 



for DHL 



(4.14) 



Pf r(7r)_ 

The physical meaning of the topological invariant in 
class A is quantized pumping of charge as a response to 
magnetic flux. In the quantum spin Hall insulator in 
class All the invariant can be interpreted either as time- 
reversal polarization pumping 26 , or as pumping of spin 
which is quantized along an unknown axis. 27,28 In the 
superconducting classes C, D, and Dili it is an analogous 
thermal or gravitational response. 29,30 



C. Topological invariant in 3D 

Turning now to 3D, we need to consider topological 
invariants of 2D Hamiltonians. The symmetry class with 
the simplest expression for the topological invariant in 
terms of the scattering matrix in 3D is AIL The 2D topo- 
logical invariant of a system in class Dili (into which All 
transforms upon dimensional reduction) is a product 25 of 
the topological invariants (4.6) of ID Hamiltonians ob- 
tained by setting one of the momenta to or 7r, 

Q[H{k u k 2 )] = Q[H(h,0)}Q[H(h,7r)], (4.15) 

with Q[H(ki)] given by Eq. (4.6). Substituting Eq. (3.4) 
into this expression we obtain 

. Pf[[/ r r(7r,0)] Vdetr(0,0) 



Pf[f/ r r(0,0)] ^detr(^,0) 



Pf [U r r{ir, tt)] ^dctr(0,7r) 
Pf[C/ T r(0,7r)] ^detr(7r,7r)' 



for AIL (4.16) 



Direct evaluation of the Hamiltonian topological in- 
variant in 2D in classes with nontrivial Chern number 
(A, C, D), and in class All is hard because of the need 
to fix the gauge throughout the Brillouin zone 12,26 . It is 
usually more efficient to use a method which relies on the 
real space structure of H evaluated in a single point in 
momentum space. 4,5,31,32 These methods using the Bott 
index or a similar expression for the topological invariant 
require the so-called band-projected position operators: 
xp = P cxp(27ria;)P and yp = P exp(2iriy)P. Here P is 
the projector on the states of the Hamiltonian with neg- 
ative energies, and a; and y are the coordinate operators 
in the unit cell of the system. In order to evaluate these 
operators in our case we note that the eigenvalues of the 
effective Hamiltonian in the symmetry classes of interest 
approach ±1 when the original system becomes localized. 
In that case P = (1 — r)/2 [with r = r(0, 0)], and we can 
avoid the need to calculate the projector explicitly if we 
approximate xp and yp by 

x P « (1 + r)/2 + (1 - r)e 27rrx {l - r)/4, (4.17) 
yp w (1 + r )/2 + (1 - r)e 27riv (l - r)/4. (4.18) 
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Using the 2D Hamiltonian expressions from Ref. 33 we 
arrive at a scattering formula for the 3D topological in- 
variant, 



Q 



1 

2^ 



Imtrlog[x P y P xpy P }, for AIII, CI, DHL 



(4.19) 

The symmetry class CII in 3D transforms upon dimen- 
sional reduction to class All in 2D. The expressions for 
the Pfafhan-Bott index required to calculate the topo- 
logical invariant for a 2D Hamiltonian in class All are 
quite involved. We do not give them here, but refer the 
interested reader to Eqs. (7), (9), and (10) of Ref. 4. 




D. Weak invariants 

All of the algorithms described above apply directly to 
the weak topological invariants. 10,34,35 In order to eval- 
uate a weak invariant one just needs to eliminate one of 
the dimensions by setting the momentum along that di- 
mension to either or 7r, and to evaluate the appropriate 
topological invariant for the resulting lower dimensional 
system. The only caveat is that since weak topologi- 
cal indices do not survive doubling of the unit cell, the 
thickness of the system in the transverse direction should 
be equal to the minimal unit cell. In the same fashion 
(eliminating one momentum or more) one can calculate 
the presence of surface states 36 in chiral superconductors 
and Fermi arcs 37 in 3D systems. 



FIG. 4: The value of the chemical potential \i c where the 
ensemble averaged topological invariant equals to 0.5, as a 
function of system size L. Red: topological invariant defined 
in terms of the scattering matrix, from Eq. (4.12). Green: 
topological invariant obtained from the Hamiltonian expres- 
sion of Ref. 4. Lines represent fits as described in the text. 



In addition to tight-binding models, our method ap- 
plies very naturally to various network models, 43-45 
which are favorite models for the phase transitions. 
Hamiltonian-based approaches are not applicable to the 
network models, since those only have a scattering ma- 
trix, and no lattice Hamiltonian. We have checked 
that calculating a topological invariant of the Chalker- 
Coddington network model of size 1000 x 1000 only takes 
several minutes on modern hardware. 



V. APPLICATIONS AND PERFORMANCE 
A. Performance 

The complexity of the Hamiltonian expressions scales 
with linear system size L as L 2 in ID, and as L 3d in 
higher dimensions. In contrast, the complexity of the 
scattering matrix expressions scales proportionally to L 
in ID and to L 3d ~ 3 in higher dimensions. 38,39 All the 
subsequent operations have the same or a more favorable 
scaling. We use the algorithm of Ref. 40 to calculate the 
Pfaffian of an arbitrary skew-symmetric matrix. 

We have verified that using the scattering matrix 
method allows to efficiently calculate the topological 
invariant of a quantum Hall system and of the BHZ 
model 41 discretized on a square lattice with a size of 
1000 x 1000. This improves considerably on previously 
reported 4,42 results of up to 50 x 50 lattice sites for the 
BHZ model. 

In 3D the improvement in performance is not as large 
because the values of L that we can reach are smaller. 
Nevertheless, we have confirmed that it is possible to cal- 
culate the topological invariant of 3D systems in classes 
All and Dili using a 4-band model on a cubic lattice with 
system size 50 x 50 x 50. This is a significant improvement 
over the 12 x 12 x 12 size, reported for Hamiltonian-based 
methods. 6 



B. Finite size effects 

The expressions for the topological invariant given in 
terms of the scattering matrix in Sec. IV do not coincide 
with Q(H) very close to the transition. This is a finite 
size effect. In order to estimate the importance of finite 
size effects we have computed the shift of the transition 
point between the n — and n = 1 plateaus of a disor- 
dered quantum Hall system as a function of size. We have 
used a square lattice discretization (lattice constant a) of 
a single band tight binding model with nearest neigh- 
bor hopping t = 1. The magnetic flux per unit cell of 
the lattice was fixed at OAh/e. We used on-site disorder 
homogeneously distributed on an interval [—0.05,0.05]. 

The transition point is defined as the value of the chem- 
ical potential fi c at which the disorder-averaged topolog- 
ical invariant equals 0.5. We have compared two expres- 
sions for the topological invariant: the scattering matrix 
expression (4.12) and the Hamiltonian expression from 
Ref. 4. The results are shown in Fig. 4. We fit the data 
obtained via the scattering matrix approach to the func- 
tion f(L) = c\ + C2/L obtaining a value C2 ~ 0.026. In 
the case of the expression of Ref. 4, the finite size effect 
arc best fit to the function g(L) = c' v +eJ> shi(e : jZ. +C4) / L, 
with c' 2 ~ 0.116. We conclude that the finite size effects 
of our algorithm are significantly lower. 
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C. Applications 

In 2D wc illustrate our approach by applying it to net- 
work models in classes A, All, and DHL In class A wc 
use the Chalker-Coddington network model. 43 In classes 
All and Dili we have used the quantum spin Hall net- 
work model of Ref. 45. In class Dili we have set the link 
phases to zero in order to ensure particle-hole symme- 
try. In each of these cases the parameter which tunes 
through the transition is the angle a related to reflection 
probability at a node of the network by R = cos 2 a. 

Our results are summarized in Fig. 5. Top panels show 
the evolution of zeros and poles of detr(z) across the 
phase transition — the 'fingerprint' of a topological phase 
transition. There are no symmetry constraints on this 
fingerprint in class A. The time-reversal symmetry en- 
sures that for every zero or pole at zq there is another one 
at 1/zq. The particle- hole symmetry translates into the 
mirror symmetry with respect to the real axis: for every 
zero or pole at zq there is one at Zq. The bottom panels 
show the behavior of the topological invariant and of the 
conductance G — titH, with t the transmission matrix 
through the system. The simulations were performed 
on systems of size 300 x 300 in each of the symmetry 
classes and averaged over 1000 samples. The presence 
of plateaus around zero in the curves for the topological 
invariant coincides with the presence of a metallic phase 
in the phase diagram of symmetry classes All and DHL 

Although wc introduced the topological invariant 
through transport properties, it does not always have the 
same features as the conductance. The topological in- 
variant characterizes winding of scattering modes in the 
transverse direction. Accordingly, in a system with a 
large ratio of width W to the length L, the width of the 
transition of the topological invariant is reduced. The 
width of the peak in the conductance, on the contrary, 
is reduced if W/L becomes small. This is in agreement 
with what we observe in numerical simulations. We have 
calculated the topological invariant and conductance av- 
eraged over 1000 disorder realizations in the Chalker- 
Coddington network model in systems with W = 300 
and L = 60 and vice versa. The results are shown in 
Fig. 6 and they agree with our expectations. 

We have also studied a 3D topological system in class 
All on a cubic lattice. We have used a simplified version 
of the Hamiltonian of Ref. 47: 



H = 



M(k) 


vk z 





vk z 


-M(k) 







vk + 


M(k) 


, vk+ 





-vk z 




— A* (5-1) 



discretized on cubic lattice with lattice constant a, where 
k± = k x ± iky, and Ai(k) = AI — ak 2 . The Hamiltonian 
parameters were chosen to be a = a 2 , v = a. We chose 
fx = fiQ + Sfj, with (iq = 0.4, and 5[i being a random mi- 
correlated variable uniformly distributed in the interval 
[—2,2]. The topological invariant defined by Eq. (4.16) 
as well as the longitudinal conductance for a 20 x 20 x 20 



system averaged over 100 disorder realizations are shown 
in Fig. 7 as a function of M. We observe that, anal- 
ogously to the two-dimensional case, the presence of a 
metallic phase is accompanied by a plateau in the topo- 
logical charge at a value of zero. 



VI. CONCLUSION 

In conclusion, we have introduced a procedure of di- 
mensional reduction which relates a scattering matrix of 
a c?-dimcnsional system to a Hamiltonian in d — 1 di- 
mensions with a different symmetry class, but with the 
same topological invariant as the original system. When 
applied repeatedly this dimensional reduction procedure 
serves as an alternative derivation of the Bott periodicity 
of topological insulators and superconductors. 

Since our approach uses only Fermi surface properties 
it is much more efficient than existing alternatives which 
require the analysis of the full spectrum. We have de- 
scribed how to implement our method efficiently in all 
the symmetry classes in 1-3 dimensions. We have veri- 
fied that it allows to analyze much larger systems than 
previously possible. 

This paper focused on the description of the method 
and we only touched on a few applications at the end. 
More applications can be envisaged and we believe that 
the scattering approach will lead to the discovery of new 
observable physics at topological phase transitions. 



Acknowledgments 

This research was supported by the Dutch Science 
Foundation NWO/FOM. We thank B. Beri, L. Fu, 
and J. Tworzydlo for useful discussions. We are es- 
pecially grateful to M. Wimmer for explaining the ef- 
ficient method to calculate transport properties and to 
C. W. J. Bcenakkcr for expert advice. 



Appendix A: Introduction to discrete symmetries 

Here we define the three core discrete symmetries, and 
the corresponding symmetry constraints on the Hamil- 
tonians and on the scattering matrices. We also specify 
how to choose the symmetry representation we used in 
Fig. 2. 



Definitions and properties of discrete symmetries 

The discrete symmetries are defined as follows: The 
time reversal symmetry operator T is an anti-unitary op- 
erator. When it is applied to an arbitrary eigenstate ip 
of the Hamiltonian H at energy e, returns an eigenstate 
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FIG. 5: Top panel: Evolution of the poles (green dots) and the zeros (red dots) of detr(z) as a function of a parameter a 
which tunes through the topological phase transition in classes A, All, and Dili in 2D. Shown is the complex plane with the 
unit circle \z\ = 1 indicated in blue. Time-reversal symmetry in All and Dili implies that for every zero/pole at jzq there is 
additionally one at 1/jzq. In Dili, there is additional particle-hole symmetry which additionally dictates zeros/poles at Zq and 
1/zo-The phase transition happens when at least one of the zeros crosses the unit circle. This event coincides with a change of 
the topological invariant Q (green) defined by Eqs. (4.12 - 4.14) 

, as shown in the bottom panels. 
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FIG. 6: Average topological invariant Q (4.12) and longitu- 
dinal conductance G of a disordered quantum Hall sample for 
different aspect ratios as a function of the mixing angle a. 



of the Hamiltonian at the same energy: 
Htp = eip => HT ip = eTip 



(Ala) 




FIG. 7: Conductance and topological invariant (4.16) for a 
disordered 3D topological insulator in class AIL 



On the other hand, the anti-unitary particle-hole symme- 
try operator V returns an eigenstate with opposite energy 
when applied to any eigenstate of the Hamiltonian: 



Hip = eip^> HVtp = -eVijj 



(Alb) 



Chiral symmetry C also reverses energy, but unlike the 
other two has a unitary operator. 

All three symmetries T, V, C are Z2 symmetries, so the 
symmetry operators must square to a phase factor. 

In an arbitrary basis the symmetry operators are rep- 
resented by 



T = U r JC, V 



U V K, 



c = u, 



(A2) 



11 



with K denoting the complex conjugation, and Uq-, U-p, 
and Uq unitary matrices. Since exp(i</>) = T 2 = Uj-U^-, 
and its determinant dct Uq-U^- is real, we obtain U-yUj- = 
±1; similarly U-pU^, = ±1. In other words, every anti- 
unitary symmetry comes in two flavors, squaring either 
to +1 or to —1. 

The symmetry constraints on the Hamiltonian 

H = THT^ 1 = U T E*\j\ (A3a) 
H = -VHV- 1 = -U v H*ul (A3b) 
H = -CHC- 1 = -U c HU f c . (A3c) 

follow immediately from the definition of the symmetries, 
Eq. (Al). 

Relation between discrete symmetries and 
translational invariance 

In addition to the basic properties, the discrete sym- 
metries in periodic systems are required to commute with 
the coordinate operator. So for any Bloch wave written 
as 

V>(r) = e* p? >(0), (A4) 

with r coordinate in a translationally-invariant system, 
and V'(O) the wave function inside a single unit cell, the 
action of the symmetry operators is: 

Ttp(r) = e- ,pr Tip(0), (A5a) 
Vtp(r) = e- tpr V^(0), (A5b) 
CV>(r) = e ipr C^(0). (A5c) 

Since the velocity of a Hamiltonian eigenstate at energy 
e and momentum p is v = de/ dp, time- reversal and chi- 
ral symmetries reverse the velocity, while particle-hole 
symmetry keeps the velocity invariant. 

Symmetry constraints on scattering matrix 

In order to figure out what the symmetry constraints 
on the scattering matrices are, we first review the ba- 
sic properties of the scattering matrices. Scattering ma- 
trices act in the space of asymptotic scattering states 
outside of the scattering region. This space contains 
two non-intersecting subspaces: the subspace of incoming 
modes and the subspace of outgoing modes. The incom- 
ing modes are all the plane waves with velocity in the di- 
rection of the scattering region, and the outgoing modes 
are all the plane waves with velocity pointing away from 
the scattering region. Let ipn be a basis in the space of 
incoming modes, and ip° nt a basis in the space of outgo- 
ing modes. Conventionally all the modes are normalized 
such that current operator in the basis of ip ln is the iden- 
tity matrix, and the negative identity matrix in the basis 
of ip out . 



The matrix elements of the scattering matrix S satisfy 

(H - e) ^V4 n + J2 S mn^ + ^ = 0, (A6) 

with ip loc a wave-function localized near the scattering 
region. 

As derived in the previous subsection, time-reversal 
and chiral symmetries change the velocity to its opposite, 
while particle-hole symmetry leaves the velocity invari- 
ant. This means that scattering states transform under 
the discrete symmetries in the following manner: 





=E(^)„ ro vc t 

m 


TC ut 






=E(^uc t ! 


cc ut 


=£(<? c )„ m «, 

m 




III 


-pi' 1 : 


hi 



(A7) 



The additional constraint on the type of time-reversal 
and particle-hole symmetries requires 

T 2 = V t Qt, V 2 = v vVv = QvQv- (A8) 

Applying time-reversal symmetry to Eq. (A6) and us- 
ing Eqs. (A7) we get 

(H - S) (v T i>° Ut + E 5 mnQTC + T^A = 0, (A9) 

where we have also used that H is time-reversal invariant. 
Comparing with Eq. (A6), we get 

S- 1 = V T 1 S*Q r , (A10) 

which we can be reduced to 

S = Q 1 -S T V r . (Alia) 

Similarly, the chiral and the particle-hole symmetry con- 
straints on S are: 

S = Q^V e , (Allb) 
S = Q i v S*V P . (Allc) 

Naturally, the constraints imposed by chiral and particle- 
hole symmetry only hold at zero excitation energy, since 
these symmetries anti-commute with the Hamiltonian. 
Finally, the symmetry constraints on the reflection ma- 
trix arc identical to Eqs. (All), since r is a diagonal 
sub-block of S. 
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Choice of symmetry representation and mapping 
from scattering matrix to Hamiltonian symmetries 

The choice of symmetry representation is fully spec- 
ified by choice of unitary matrices Uq, Vq, and Qo 
(O = T,V, or C). The symmetry representations used 
in the main text were chosen to make the mapping from 
the reflection matrix to an effective Hamiltonian most 
straightforward. In order to reach this aim, we always 
choose Vo = ±Qo for each of the three symmetries. 
Whenever chiral symmetry is present, we use 

U c =t z , V c = Qc = l. (A12) 

Our choices of V and Q for the other two symmetries 
with 

V r = V v = V, Q r = Qv = Q (A13) 

depend on the specific symmetry class. When V 2 = — 1 
(symmetry classes C, CI, CII) we choose V = a y , and 
we choose V = 1 in symmetry classes D, Dili, and BDI, 
where V 2 = 1. The relative sign between V and Q follows 
from Eq. (A8): 

q = v 2 T 2 V. (A14) 

In the remaining two classes AI and All we choose V = 1 , 
and Q = T 2 ■ Symmetry representations of the effective 
Hamiltonians follow immediately from Eqs. (3.4). 

Finally we show how symmetry operators change upon 
creating an effective Hamiltonian from a reflection ma- 
trix. The effective Hamiltonian created from a reflection 
matrix with chiral symmetry satisfies 

H = V 2 T 2 VH*V, (A15) 

so that the resulting symmetry of the effective Hamilto- 
nian is particle- hole if V 2 T 2 — — 1, and time- reversal if 
V 2 T 2 = 1. The symmetry operator of this symmetry 
squares to VV* . If a reflection matrix has only time- 
reversal symmetry, then the time-reversal and particle- 
hole symmetry constraints on the effective Hamiltonian 
have the form 

H = ±V®t x H*V®t x , (Al6a) 
H = tV® r y H*V® T y , (A16b) 



where the sign is determined by the choice of representa- 
tion of the symmetry of r. Hence the symmetries of the 
effective Hamiltonian satisfy T 2 = —V 2 . Finally, for an 
effective Hamiltonian constructed from a reflection ma- 
trix with only particle-hole symmetry, the resulting sym- 
metry constraints on the Hamiltonian are 

H = ±V®t H*V ® t , (Al7a) 
H = TV® t z H*V® t z , (A17b) 

so that both symmetry operators square to the same 
value. 

Appendix B: Calculation of the number of poles of 
detr(z) inside the unit circle 

We prove that the equation 

det[A-Zi(k)]=0 (Bl) 

has Ni solutions with z = e lk inside the unit circle, and 
Ni solutions outside of the unit circle as long as A' A 
only has eigenvalues less than one, which is generically 
the case since A is a sub-block of a unitary matrix S. Let 
us assume that ip is an eigenvector of the corresponding 
eigenvalue problem: 

Ail> = Zt(k)il>, (B2) 

with an eigenvalue with \z\ = 1. In this case 
Zl(k)Zi(k) = 1. We come to a contradiction by con- 
sidering the following inequality: 

ipty > i^A^Ail) = ^ Z\{k)Zi(k)i) = ipty (B3) 

So we conclude that there are no solutions of det[A — 
Z\(k)] = with z on the unit circle. Next, we observe 
that for A = there are exactly Ni of dctLA— Z\{k)\ = 
with z = and N\ solutions with z = oo. Since these 
solutions never cross the unit circle when A is smoothly 
deformed, we come to the statement we set to prove. 
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